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Abstract — Cellular systems are becoming more heterogeneous 
with the introduction of low power nodes including femtocells, 
relays, and distributed antennas. Unfortunately, the resulting 
interference environment is also becoming more complicated, 
making evaluation of different communication strategies chal- 
lenging in both analysis and simulation. This paper suggests 
a simplified interference model for heterogeneous networks. 
Leveraging recent applications of stochastic geometry to analyze 
cellular systems, this paper proposes to analyze downlink perfor- 
mance in a fixed-size cell in what is called the hybrid approach. 
The interference field consists of a contribution from out-of-cell 
interferers modeled as a superposition of marked Poisson point 
processes outside a guard region and a cross-tier contribution 
modeled as a marked Poisson point process. The guard region is 
used to avoid association issues that would occur if an interferer 
were located at the cell edge. Bounding the interference power as 
a function of distance from the cell center, the total interference 
is characterized through its Laplace transform. An equivalent 
marked process is proposed for the out-of-cell interference under 
additional assumptions. To facilitate simplified calculations, the 
interference distribution is approximated using the Gamma 
distribution with second order moment matching. The Gamma 
approximation simplifies calculation of the success probability 
and ergodic rate, incorporates small-scale and large-scale fading, 
and works with co-tier and cross-tier interference. Simulations 
show that the guard region can be tuned to give a reasonable 
match with performance in a hexagonal grid, and that the 
Gamma approximation provides a good fit over most of the cell. 



I. Introduction 

Cellular network deployment is taking on a massively 
heterogeneous character as a variety of infrastructure is being 
deployed including macro, pico, and femto base stations 1 1 1, as 
well as fixed relay stations |2| and distributed antennas |3|. A 
major challenge in deploying heterogeneous cellular networks 
is managing interference. The problem is compounded by 
the observation that the deployment of much of the small 
cell infrastructure will be demand based and more irregular 
than traditional infrastructure deployments ||4|. Further, as 
urban areas are built out, even macro and micro base station 
infrastructure is becoming less like points on a hexagonal 
lattice and more random |[5|. Consequently, the aggregate 
interference environment is more complicated to model and 
evaluating the performance of communication techniques in 
the presence of heterogeneous interference is challenging. 

Stochastic geometry can be used to develop a mathematical 
framework for dealing with interference (ISll-llQl. With stochas- 
tic geometry, interferer locations are distributed according to 
a point process, often the homogeneous Poisson Point Process 
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(PPP) for its tractability. PPPs have been used to model co 
channel interference from macro cellular base stations ||5 
1 10|-|fT3l, cross-tier interference from femtocells Q, |14 
1 15 1, co-channel interference in ad hoc networks ||7|~l|D' 
and as a generic source of interference fT7{-|fT9[. 
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Modeling co-channel interference from other base stations 
as performed in g, |[T0), |[12), ||T3), ||20) is a good starting 
point for developing insights into heterogeneous network inter- 
ference. In 1 10 1 PPPs are used to model various components 
of a telecommunications network including subscriber loca- 
tions, base station locations, as well as network infrastructure 
leveraging results on Voronoi tessellation. In [12] , [20] a 
cellular system with PPP distribution of base stations, called a 
shotgun cellular system, is shown to lower bound a hexagonal 
cellular system in terms of certain performance metrics and 
to be a good approximation in the presence of shadow fading. 
In [31, a comprehensive analysis of a cellular system with 
a PPP is presented, where key system performance metrics 
like coverage probability and average rate are computed by 
averaging over all deployment scenarios. The approach in Q 
is quite powerful, leading to closed form solutions for some 
special choices of parameters. In p3) the outage and the han- 
dover probabilities in a PPP cellular network are derived for 
both random general slow fading and Rayleigh fading, where 
mobiles are attached to the base station that provides the best 
mean signal power. Recent work has considered extensions 
of 1 5), p2) to multi-tier networks. For example \2 1 J extends 
| [T2 to provide some results on the signal-to-interference 
ratio in multi-tier networks while 1 22 1 extends fS] to provide 
remarkable simple expressions for the success probability in 
multi-tier networks. From our perspective, the main drawbacks 
of the approach in |5|, |10|, |12|, |13), |j20| for application 
to existing systems is that performance is characterized for 
an entire system not for a given cell. As cellular networks 
are already built out, a cellular provider will often want to 
know the performance they achieve in some given cells by 
adding additional infrastructure (and thus interference) in the 
rest of the network. It is also challenging to incorporate more 
complex kinds of heterogeneous infrastructure like fixed relays 
or distributed antennas into the signal and interference models. 

Models for co-channel interference from Poisson Point Pro- 
cesses in general wireless settings has been considered in prior 
work |(l6l-|[T9l, (23). References |[17)-|[19;| derive models for 
the complex noise distribution. These models assume that the 
noise changes on a sample-by-sample basis and are suitable 
for deriving optimum signal processing algorithms for example 
optimum detectors. Other work, see e.g. y_6J, ||23| and the 
references therein, provide models for the noise power distri- 
bution. Most system-level analysis work fTO) , |12|, p3) , 



pO) focuses on the noise power distribution. They assume 
that the interference distribution conditioned on the noise 
power is Gaussian and thus characterize performance based 
on quantities like the signal-to-interference ratio or the signal- 
to-interference-plus-noise ratio, which are not appropriate for 
other noise distributions [17J-[19J. This is reasonable because 
cellular systems tend to have structured transmissions in time 
and frequency, the interfering sets are likely to be constant 
over a coding block. Consequently we focus on characterizing 
and employing the noise power distribution in our work. 

In this paper we propose a simplified interference model for 
heterogeneous networks called the hybrid approach. The key 
idea is to evaluate performance in a fixed-size cell considering 
co-channel interference outside of a guard region and cross- 
tier interference from within the cell. Different sources of 
interference are modeled as marked Poisson point processes 
where the mark distribution includes a contribution from 
small-scale fading and large-scale fading. The guard region is 
used to avoid association issues that would occur if co-channel 
interferer were located at the cell edge. Based on simulations, 
an optimum guard region size of two-thirds the cell radius is 
suggested to achieve good match with hexagonal interferer 
locations. For co-channel interference, an upper bound on 
the interference power distribution is characterized through its 
Laplace transform as a function of distance from the cell center 
for the case of downlink transmission, after lower bounding the 
size of the exclusion region. When there are multiple types of 
co-channel interference, the distribution power can be further 
upper bounded by a single marked process that is a function 
of the mark parameters, densities, and minimum guard radius 
from each interference process. One of the key advantages of 
the proposed model is that no explicit user association scheme 
needs to be specified. Most prior work on heterogeneous 
cellular networks analysis using stochastic geometry derives 
the SINR distribution assuming a specific user assignment 
scheme. In | |24| , the downlink SINR distribution from the 
serving BS to an arbitrarily-located user is derived when the 
candidate serving BS at each tier is the nearest to the user. 
In ||25l the user is attached to the BS providing the strongest 
average received signal, whereas in p2| , p6j , user association 
based on maximum SINR connectivity is considered. In the 
hybrid model, users can be arbitrarily distributed in the cell 
of interest without requiring any assumption on the user 
association. This allows us to analyze the cell performance 
and calculate per-user performance metrics under various user 
distributions. 

To facilitate simplified calculations, we use approximations 
based on the Gamma distribution. We model the small-scale 
fading component using the Gamma distribution as this incor- 
porates single user and multiple user signaling strategies with 
Rayleigh fading as special cases. We approximate the product 
of the Gamma distributed small-scale fading component and 
the log-normal large-scale fading component as a Gamma 
distributed random variable, based on prior work p7) . We then 
approximate the interference distribution using the Gamma 
distribution with second order moment matching; the moments 
are finite because we use a nonsingular path-loss model. Using 
the Gamma approximation, we present simplified expressions 



for the success probability and the average rate, assuming 
that the system is interference limited, i.e. our calculations 
are based on the signal-to-interference ratio. Simulations show 
that the Gamma approximation provides a good fit over most 
of the cell, with larger error towards the cell edge. 

Our mathematical approach follows the approach of ||5| 
and leverages basic results on PPPs ||8], ||9|, |23|. Compared 
with fSl, we consider a fixed-size cell and do not evaluate 
performance in a Voronoi tessellation of randomly distributed 
interferers. We also do not calculate system- wide performance 
measures, rather we focus on performance for a fixed cell of 
interest. We demonstrate how to employ the proposed fixed 
cell analysis in a heterogeneous network consisting of mixtures 
of the different kinds of infrastructure. The advantage of our 
approach is that more complex types of communication and 
network topologies can be analyzed in a given target cell and 
key insights on the performance of advanced transmission 
techniques with heterogeneous interference using stochastic 
geometry can be obtained, at the expense of generality. Com- 
pared with our preliminary results reported in f28l, in this 
paper we deal with more general fading distributions and use 
the Gamma approximation to simplify the calculation of the 
success probability and the ergodic rate. Guard zones have 



been used in other work on interference models e.g. 1 19 1, p9) . 
In | |29| they are used to evaluate the transmission capacity in 
contention-based ad hoc networks not cellular networks while 
in 1 19] they are considered in the derivation of the baseband 
noise distribution but not the noise power. 

II. Mathematical Preliminaries 

In this paper we make frequent use of Gamma random 
variables. We use them to model the small-scale fading distri- 
bution, to approximate the product of the small-scale and log- 
normal fading distribution, and to approximate the interference 
power. In this section we provide some background on Gamma 
random variables to make the paper more accessible. Proofs 
are omitted where results easily follow or are well known in 
the literature. 

Definition 1: A Gamma random variable with finite shape 
fc > and finite scale > 0, denoted as r[fc, 6] has probability 
distribution function 



fx{x) = X 



6''=r(/fc) 

for x > 0. The cumulative distribution function is 

r(fc) 



Fx{x) 



(1) 



(2) 



where 7(a, h) = *dt is the lower incomplete gamma 

function. The first two moments and the variance are 



EX = ke EX^ = fc(l + fc)6'2 var X 



k0^. 



(3) 



The scale terminology comes from the fact that if X is r[fc, 6] 
then with scalar a > 0, aX is r[fc, ad]. 

The Gamma distribution is relevant for wireless communi- 
cation because it includes several channel models as a special 
case. For example if Z is a circularly symmetric complex 
Gaussian with Af{0,l) (variance 1/2 per dimension) then 
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X = |Zp is r[l, 1], which corresponds to Rayleigh channel 
assumption. If X has a Chi-square distribution with 2N 
degrees of freedom, used in the analysis of diversity systems. 



then X is also r[iV, 2]. If X is Nakagami distributed (|30| 
with parameters /i and uj then Y = X^ has a r[fc,6'] 
distribution with fc = /i and 6 — Lo/ji. Consequently the 
Gamma distribution can be used to represent common fading 
distributions. 

Of relevance for computing rates, the expected value of 
the log of a Gamma random variable has a simple form as 
summarized in the following Lemma. 

Lemma 2 (Expected Log of a Gamma Random Variable): 
If X is T[k,e] then ElnX = i}{k) + ln(6l) where ^{k) is 
the digamma function. 

The digamma function is implemented in many numerical 
packages or for large x its asymptotic approximation can be 
used '4){x) = \n{x) + l/x + 0{l/x'^). 

The Gamma distribution is also known as the Pearson type 
III distribution, and is one of a family of distributions used to 
model the empirical distribution functions of certain data |3ll. 
Based on the first four moments, the data can be associated 
with a preferred distribution. As opposed to optimizing over 
the choice of the distribution, we use the Gamma distribution 
exclusively because it facilitates calculations and analysis. 
We will approximate a given distribution with a Gamma 
distribution by matching the first and second order moments 
in what is commonly known as moment matching p2| . The 
two parameters of the Gamma distribution can be found by 
setting the appropriate moments equal and simplifying. We 
summarize the result in the following Lemma. 

Lemma 3 (Gamma 2""^ Order Moment Match): Consider a 
distribution with /i = 'EX, /i*^^'' = EX^, and variance = 
/i*^^' — /i^. Then the distribution T[k,9) with same first and 
second order moments has parameters 



k = 



(4) 



Note that the shape parameter is scale invariant. For exam- 
ple, the Gamma approximation of aX would have shape 
k = /^^/c^ and scale 6 — aa^ j [i. 

We will need to deal with sums of independent Gamma 
random variables with different parameters. There are various 
closed-form solutions for the resulting distribution p3), (34| 



(see also the references in the review article |35 1). In this paper 
we choose the form in | |33| as it gives a distribution that is a 
sum of scaled gamma distributions. 

Theorem 4 (Exact Sum of Gamma Random Variables): 
Suppose that are independent T[ki,9k] distributed 

random variables. Then the probability distribution function 
of y = Yi can be expressed as 



fYiy) = C 



(5) 



mill, 0i, p = Yl"=l C" = n"=l(^rnm/^i) 



where O^i 
and Ci is found by recursion from 

1 



ki 



i+l 



1 



and 7„j = J27=i - 0min/di)"^- 

Proof: See f33l. ■ 

For practical implementation, the infinite sum is truncated; 
see [33) for bounds on the error. Mathematica code for 
computing the truncated distribution is found in (36). In some 
cases, many terms need to be kept in the approximation to 
achieve an accurate result. Consequently we also pursue a 
moment-matched approximation of the sum distribution. 

Lemma 5 (Sum of Gammas 2"'' Order Moment Match): 
Suppose that {Yi} are independent Gamma distributed 
random variables with parameters ki and 9i. The Gamma 
distribution r[ky,6y] with the same first and second order 
moments has 



and 6.,, 



(7) 



Bounds on the maximum error obtained through a moment 
approximation are known |32| and may be used to estimate the 
maximum error in the cumulative distribution functions of the 
Gamma approximation. In related work fyT\ we show that the 
approximation outperforms the truncated expression in ( fTT] ) 
unless a large number of terms are kept, and therefore conclude 
that the approximation is reasonable. The approximation is 
exact in some cases, e.g. if all the shape values are identical 



then ky = h 



9rr„r, and 



(in this case 9i 

C^ = 0). 

In wireless systems, the zero-mean log-normal distribution 
is used to model large-scale fluctuations in the received signal 

(Inx))^ 

power. The distribution is fx{x) = ^ e 2<t2 where 

a is the shadow standard deviation, often given in dB where 
= 8.686(7. Typical values of cr^s are between 3dB and 



lOdB for example in 3GPP 1 38 1, 1 39 1 



Composite fading channels are composed of a contribution 
from both small-scale fading and large-scale fading. Consider 
the random variable for the product distribution HL where 
H is r[kfi,dh] and L is log-normal with parameter a. Related 
work has shown that the product distribution can be reasonably 
modeled using a Gamma distribution |27| by matching two of 
the moments. In our case we find that the first and second 
central moments gives reasonable performance. 

Lemma 6 (Gamma Approximation of Product Distribution): 
Consider the random variable PHL where P is a constant, H 
is T[kh, 9h\ and L is log-normal with variance a. The Gamma 
random variable with T[kp,Op] with the same first-order and 
second-order moments has 

(8) 



kp — 



{l/kh + 
{l + kh) 



l)e'^' - 1 



khOhe' 



72 



Proof: Note that 

EHL= EH EL = khOhe"^/^ 



and 



var HL = (1 + khjkhOle^''^ 



(9) 



(10) 



(11) 



(6) Substituting into Lemma [3] and simplifying yields ([H) and (|9|. 
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III. Downlink Network Model 

In the classic model for cellular systems in Fig. [TJa), base 
stations are located at the centers of hexagons in a hexagonal 
tessellation and interference is computed in a fixed cell from 
multiple tiers of interferers. The hexagonal model requires 
simulation of multiple tiers of interferers and makes analysis 
difficult. In the stochastic geometry model for cellular systems 
illustrated in Fig.[TJb), base stations positions follow a PPP and 
cells are derived from the Voronoi tessellation. Cell sizes are 
random and performance metrics are computed in an aggregate 
sense accounting for all base station distributions. Our pro- 
posed model, which we call the hybrid approach, is illustrated 
in Fig. [TJc). We consider a cell of fixed shape and size, 
nominally a circle with radius Re- The base station locations 
outside of the fixed cell are modeled according to a PPP and 
the interference becomes a shot-noise process pO) . Since the 
edge of the cell is fixed, and is not determined by interference 
levels, a guard region of radius Rg from the cell edge is 
imposed around the fixed cell in which no other transmitters 
can occupy. Performance is evaluated inside the fixed cell 
accounting for interference sources outside the guard region. 
In general, Rg is a parameter of the model; our simulations 
show that Rg — 2Rc/3 provides good performance relative to 
a hexagonal tessellation. 




(a) Fixed Geometry (b) Stochastic Geometry (c) Hybrid Approach 



Fig. L Models for cellular communication, (a) The common fixed geometry 
model with hexagonal cells and multiple tiers of interference, (b) A stochastic 
geometric model where all base stations are distributed according to some 2D 
random process, (c) The proposed hybrid approach where there is a fixed cell 
of a fixed size surround by base stations distributed according to some 2D 
random process, possibly with an exclusion region around the cell. 

To make the calculations concrete, in this paper we focus 
on the downlink. We consider the received signal power at 
distance r from the transmitter presumably located at the cell 
center. Extension to distributed antennas for example requires 
a more complex signal model. Let the received signal power 
be 

S{r) = P{r)LH, (12) 

where P{r) = Ps/l{r) is the distant-dependent average 
received signal power with Pg being the transmit power, L a 
random variable corresponding to the large-scale fading power 
usually log-normal, and is a random variable corresponding 
to the small-scale fading power. We assume that a non-singular 
path-loss model is used with 

^(r) = Cmax(do,r-)", (13) 

where n > 2 is the path-loss exponent, do > is the 
reference distance, and C > is a constant. In addition 



to resulting in finite interference moments, ( [T3] l takes into 
account the realistic RF design constraints on the maximum 
received power and is often assumed in a standards based 
channel model like 3GPP LTE Advanced i38J. 

We consider the power due to small-scale fading as oc- 
curring after processing at the receiver, e.g. diversity com- 
bining. Furthermore, we assume that the transmit strategy is 
independent of the strategies in the interfering cells and no 
adaptive power control is employed. We model the small- 
scale fading power H as a Gamma distributed random variable 
with r[kh,Oii]. The Gamma distribution allows us to model 
several different transmission techniques. For example, in con- 
ventional Rayleigh fading H is r[l, 1] while Rayleigh fading 
with A^i receive antennas and maximum ratio combining H is 
r[A'i., 1]. Maximum ratio transmission with Rayleigh fading, 
Nt transmit antennas, and = 1 receive antennas would 
have H is r[Nt, 1]. With multiuser MIMO (MU MIMO) with 
A't transmit antennas, U < active users, and zero-forcing 
precoding then H is T[Nt-U+l,l/U] where the l/U follows 
from splitting the power among different users. 

We model the large-scale fading contribution L as log- 
normal distributed with parameter a^. The large-scale fading 
accounts for effects of shadowing due to large objects in the 
environment and essentially accounts for error in the fit of 
the log-distance path-loss model. Because of the difficulty 
in dealing with composite distributions, we approximate the 
random variable HL with another Gamma random variable 
with the same first and second order moments as described in 
Lemma |6] with distribution r[kp,6p] where kp is found in ^ 
and dp in (|9]l. Note that from the scaling properties of Gamma 
random variables, with this approximation, the received signal 
power S{r) is equivalently r[kp, P{r)6p]. 

IV. Homogeneous Network Interference 

In this section we develop a model for homogeneous 
network interference, i.e. interference comes from a single 
kind of interferer such as a macrocell. We model the inter- 
ferer locations according to a marked point process, where 
the marks correspond to the channel between the interferers 
and the target receiver. Specifically we consider a PPP $ 
with density A and marks modeled similarly to the signal 
power distribution. For interferer k, Lk follows a log-normal 
distribution with a while the small-scale fading distribution Gk 
is V[kg^9g]. The transmit power of the interferer is denoted 
by Pg. The transmit power is fixed for all interferers and 
the actions of each interferer are independent of the other 
interferers. Note that the parameters of the interference mark 
distribution are usually different than the signal channel even 
when the same transmission scheme is employed in all cells. 
For example, with single antenna transmission and conven- 
tional Rayleigh fading Gk is r[l, 1] for any N,. at the receiver, 
while with multiuser MIMO zero-forcing beamforming with 
multiple transmit antennas serving U < Nt users Gk is 

r[c/,i]. 

We aim at computing the interference power received at 
a distance < ?' < i? away from the center of the cell as 
illustrated in Fig. |2] Unfortunately, the exclusion distance to 
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Fig. 2. Interference power calculation at a point r away from the transmitter. 
An upper bound is considered as the interference is not excluded uniformly 
around the receiver. 



the nearest interferer is asymmetric, since the distance to the 
closest edge of the circle is Rc + Rg — r while the distance to 
the furthest edge is Rc+Rg+r. To avoid the dependence on the 
location of the interference field we pursue an upper bound 
on the interference power, as illustrated Fig. |2] Specifically, 
we consider the interference contribution in a ball of radius 
Rc+Rg—r around the received signal denoted B{Rc+Rg—r). 
This results in an upper bound on the aggregate interference 
due to the fact that less area is excluded from the calculation. 
Then the receiver is at the center of the reduced interference 
region, and because it is a PPP thanks to Slivnyak's theorem 
Q can be assumed to be at the origin. Under this assumption 
we write the received power as 



/(r) 



E 

fe6*\B(i?.<:+i?.g-l-) 



GkLkPg 



(14) 



We use the same non-singular path-loss model in ( [T3| l as 
for the transmit signal power. Other non-singular models and 
tradeoffs between different models are discussed in |8J, | ,23J . 
As long as there is a guard region with i?g > a singular path- 
loss model could be employed without changing the results. In 
our calculations we assume that Rc + Rg — r > to simplify 
the exposition. 

We will characterize the distribution of the random variable 
I[r) in terms of its Laplace transform, along the lines of fs). 
When the interference marks GkL^ follow an arbitrary but 
identical distribution Vfc the Laplace transform is given by 



£,(,)(s)=E,(,)[e-^^W] 



Ed 



fce«\B(fic + Rg 



yk(iiS>\B(R„+R.g-r) 

E$ I W Eg,L (e^'W 

fee*\B(i?,,+fl,g-r) ^ 

(15) 



(a) 



--F,i.Rc + Hg->-.c(s) 



(16) 



where (a) follows from the i.i.d. distribution of G^Lk and its 
further independence from the point process $, [h] follows 
assuming that Rc + Rg — r > do and using the probability 
generating functional of the PPP fE\, Z — GkL^Pg, and 



Fn. 



x,C 



n 



(17) 



In some cases the expectation can be further 
evaluated, e.g. when Z ^ ^[k,0] then F„.2;,c(s) = 



2+kn 



FAk,k 



A + k- 



-6'2/"6(fc+ where B{x,y) is the Beta Euler 
function and 2P1 is the Gauss hypergeometric function. 



V. Heterogeneous Network Interference 

In this section we extend the proposed interference model 
to the case of a heterogeneous network, where the macrocel- 
lular network is underlaid with other kinds of infrastructure 
including low-power nodes like small-cells, femtocells, fixed 
relays, and distributed antennas. The key idea is to model 
the interference from each type of infrastructure separately 
using a marked PPP described by the model in Section |IV] 
and then develop an equivalent marked PPP for the sum 
interference term to facilitate analysis. We consider two dif- 
ferent interference scenarios as illustrated in Fig. |3] The 
scenario in Fig. [3ja) corresponds to multiple kinds of out- 
of-cell interference. For example, there might be interference 
from both high power base stations and low power distributed 
antenna system transmission points. The guard regions associ- 
ated with different processes may be different. The scenario in 
Fig. [3jb) corresponds to cross-tier interference from a second 
tier of low-power nodes, e.g. small-cells or femtocells or 
uncoordinated transmission points. The main difference in this 
case is that there is interference within the cell so effectively 
the interference contribution is a constant since it is spatially 
invariant. There may still be a guard radius in this case but it 
is likely to be small. 




(a) (b) 



Fig. 3. Interference in a heterogeneous network from multiple kinds of 
infrastructure, (a) Heterogeneous out-of-cell interference, (b) Homogeneous 
interference and cross-tier interference from low power nodes. 
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A. Heterogeneous Out-of-Cell Interference 

Suppose that there are AI different kinds of infrastructure 
deployed in a homogenous way through the network. A more 
complex model would take non-uniformity or clustering into 
account, but this is beyond the scope of this work. Each 
interference source can be modeled as a marked PPP with 
marks corresponding to the composite fading distribution 
distributed as r[fc,„, 9,n] where m = 1, . . . , M. The path-loss 
function is the same for each process. The transmitting power 
is given by Pm, the transmitter density by Am, and the guard 
radius by i?g™\ Consequently, each process is parameterized 
by (fcm, 6'm, Pm, Am, This gives an equation for the 

total interference as 



Ju 



E 



™-i fee$„AS(flc+-Rg'"'-'-) 



^{rn) J- (m) p 
^fc ^" 



(18) 



The difference in guard regions adds a complication to 
developing an equivalent marked PPP. Consequently we up- 
perbound the interference level as described in the following 
lemma. 

Lemma 7: Consider the sum of m PPPs with parame- 
ters {kjyi^ Om-j 

Pm,Am,4™^ Then 



M 



JUr) < E E 

m=l fee$„\S(_Re+-R„ 



^(m) j(m) p 
'-'k ^k 



J{r) 



(19) 



chosen point (typical point) of a; G $s, then the probability 
for the event that x belongs to $m \ B{Rc + i?min — i"). 



i.e. P(xe $m\S(P, + i?,„i„-r)) 

.(e) 



Hence, 



the random variable X^, (interference mark) of any 
point X e $s is equal to G'^'^^L^'^^Pm with probability 
P (a; e $m \ B{R, + i?„,i„ - r)) = Am/ A„. ■ 

The marks of the cumulative heterogeneous interference has 
a finite mixture distribution fx{x) given by 



Z^m=l 



(21) 



where /^""^(z) denotes the PDF of = G^"''^ L^'^^ Pm- 

From the independence of the M PPPs and using similar 
arguments as for deriving ( fTS] ), the Laplace transform of the 
heterogeneous out-of-cell interference is 



M 



(22) 



m— 1 



with £j(m)/^(s) the Laplace transform of j''™'\r) defined as 
in Lemma ^ with appropriate substitution for Pm, Am, and 

Expanding the integrals in ( |22| i along the lines of the 
calculations that lead to ([16]), combining terms, and using the 
equivalent mark distribution from Proposition [8] leads to the 
following result. 

Corollary 9: The Laplace transform of the cumulative in- 
terference from all M tiers is given by 



where Pmin = min Pf,™-* and X < Y denotes that X is 

l<m<J\/ " 

Stochastically dominated by Y. 

Proof: Follows directly by noting that 



(m) 



vol i^B{Rc + R 
and that $m \B{Rc + R. 



(m) 
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> vol {B{R, + R^in-r)) 
-r) C $m\S(Pc + i?mi„-r-). 



The superposition of M different homogeneous PPPs can 
be treated as a single equivalent PPP by modifying appropri- 
ately the mark distribution. This many simplify analysis and 
simulations in some cases. 

Proposition 8: The cumulative interference received from 
all M interfering networks 



M 



^(m) j-{m) p 
^k ^k 



E E 

m=l fce$,„\B(fle+flmin-r') ^("^k ) 



.(e) 



\ ^ x^ 



(20) 



where = | |J $m \ S(Pc + i?mi„ - r-)| and xj^'^ is a 
random variable that is equal to G''"''L'^'"'Pm with probability 



^rn/ J2m=l ^ra for TO = 1, 



Proof: The union of M PPPs, i.e. $5 ^ 
{ U ^m\B{Rc + Rimn- r)\, is a PPP with intensity 

measure Ae = X)m=i since the superposition of M 
independent PPPs is also a PPP |8J. Consider now a randomly 



Cj{s) = e^eTr(i?c+_R.„,„-r)^-?I^iiii^^i^F„,H, + H„„„-.,c(s) 

(23) 

where Ag = I]m=i Am- 
using Corollary [9] we can replace calculations involving the 
heterogeneous interference term with an equivalent marked 
interference process. Consequently the homogeneous inter- 
ference results can be applied but with a different mark 
distribution and a total density term. 

B. Cross-Tier Interference 

Certain kinds of infrastructure, like femtocells, are not 
associated with the macro base station. Rather they form 
another tier of nodes and create cross-tier interference. We 
propose to use the same stochastic geometry framework to 
model the interference with a main difference in how the 
notion of guard zone is defined. Let Bcross denote an exclusion 
region around the desired receiver with radius Rcross- This 
region might be quite small, just a few meters radius, and is 
designed to avoid the case where the low-power node is co- 
located with the target receiver Let us suppose that the cross- 
tier interference is associated with a PPP given by Across, 
the transmitting power is given by Pcross, and the transmitter 
density by Across- The received interference power is given by 

GkL^Pcross 



P 



E 



t{Rk 



(24) 
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which does not depend on r due to the shift invariance property 
of the PPP. Consequently cross-tier interference creates a 
constant interference that is independent of the mobile receiver 
location. 

VI. Approximating the Interference Distribution 
Using the Gamma Distribution 

The expressions in Section |IV] and Section |V] provide a 
characterization of the distribution of the interference terms in 
the proposed hybrid approach. Unfortunately, the distributions 
are characterized in terms of the Laplace transform and further 
simplifications exist only for special choices of the parameters 
and mark distribution. Consequently, in this section, we pursue 
an approximation of the distribution of the interference using 
the Gamma distribution. First we review the calculations for 
the Gamma approximation of the interference terms with a 
single and multiple interferers. Then we use these calculations 
to parameterize the Gamma approximation of the interference. 

A. Moments of the Interference Process 

The moments for interference distribution can be computed 
along the Unes of (Sj Equations 2.19 and 2.21]. We only 
need the first two moments for the Gamma approximation. Let 
Z = GkLkPg denote the random variable corresponding to the 
mark distribution for the case of a homogeneous interference 
source. Then using results from fSl 



E/(r) = (EZ) 2tt\ 



E/(r)2 = (E/(r))2 



Ra+ii„ 



w 



dw 



POO 

{EZ^) 2ttX i 



l{w) 



(25) 



Tdw. 



(26) 

The interference is a type of shot noise. For the nonsingular 
model of ([T3|, the value is finite for any r. Now we assume that 
Rc + Rg — r > d[) to simplify the exposition. We summarize 
the results in the following proposition. 

Proposition 10 (Moments of Homogeneous Interference): 
The mean and variance of the interference in the case of 
homogenous interference are 

[R, + Rg- r)^-" 
C(n - 2) 



EI{r) = (EZ) 27rA- 



(27) 



var I{r) ^ (EZ^) 2-k\ 



(i?c + Rg 



\2-1n 



(28) 



2C2(n- 1) 

In the case of heterogeneous interference, we propose 
to use the moments of the sum to approximate using the 
equivalent mark distribution. We summarize the key results 
in the following proposition. 

Proposition 11 (Moments of Heterogeneous Interference): 
The mean and variance of the interference in the case of 
heterogeneous interference are 



Rn 



-rf 



m— 1 ^ ' 

M , 



(29) 



Ug 



2C2(n- 1) 



where = g'-^^ l'^^ Pm- 

The calculations can be further simplified using the minimum 
radius approximation but we defer their presentation due to 
lack of space. 

The approach in ( [TT| is quite general. An additional term 
could be included accounting for the cross-tier interference 
Icross in a straightforward way. Further modifications are 
possible. For example, having different numbers of active 
users in an interfering cell with multiuser MIMO could be 
included by having a sum of marked point processes each 
with a different fading distribution corresponding to different 
choices of active users. Note that we can incorporate cross-tier 
interference Icross into Proposition [TT|by repl acing Rq^Rq — t 
with Rcross and selecting the other parameters as appropriate. 

B. Gamma Approximation of the Interference 

While the Laplace transform of the interference power 
completely characterizes the distribution, it is often chal- 
lenging to provide simple or closed-form solutions for a 
variety of scenarios of interest. Consequently we pursue an 
approximation of the interference term that will yield simpler 
expressions targeted at the proposed hybrid setting using the 
Gamma distribution. 

Proposition 12 (Homogeneous Interference): The r[A:i,0j] 
random variable with the same mean and variance as I{r) 
has 

.{EZf 2{n - 1) 



ki = 2ti\{Rc 



Rg^r) 



EZ2 {n-2f 



and 



EZ^ 
EZ 



{Rc 



Rg-r) 




(31) 



(32) 



kjOi 



Proof: Using the results of Lemma p 
and var /(r) — ki9^, then substituting ([^ i and ( [28] l, gives ki 
and 9i. ■ 
We can derive a similar expression in the case of heteroge- 
neous interference, using Lemma [3] We omit the result due to 
lack of space. 

Now suppose that the mark distribution PgL^Gk is approx- 
imated by a Gamma distribution. Substituting and simplifying 
leads to the following corollary. 

Corollary 13 (Homogeneous Interf. - Gamma Fading): 
Suppose that the product fading distribution LkGk is 
T[kz;Oz\. Then the parameters of the second order Gamma 
approximation are 

-47rA(i?e 



1 



Rg~ry 



(n-l) 



and 



e, = Pge,{i + k,){Rc + Rg-ry 



{n-2Y 



1 n 



(33) 



(34) 
then 



EZ 



2Cn-V 

Proof: If Z is PgL^Gk where L^Gk is T[kz 

and EZ^ = k^P^9l{l + k^). Substituting into 
pT| ) and ( |32| l gives the result. ■ 
We also provide the result for the heterogeneous case where 
each tier of interferes has power f,„ and product fading dis- 
tribution given by G^.^'l'''" 



f, with distribution r[fcm, ^m^W 



^^^^ Substituting and simplifying gives the second corollary. 



7 



Corollary 14 (Heterogeneous Interf. - Gamma Fading): 
Suppose there are M tiers of interference, each with power 
Pm and product fading distribution given by L^™'' 
approximated as T\k„i,9mPm]. Then the parameters of the 
second order Gamma approximation are 



and = ^ 
varJ(r) EJ(r) 



(35) 



where 



M 



m— 1 

M 



C{n - 2) 



(36) 



var J (7 



m— 1 



(i?c + i?: 



(m) 



2C2(n- 1) 



(37) 

Proof: The main equation follows from Lemma |3] insert- 
ing the moments for each interference tier computed along the 
lines of Proposition [lO] ■ 

VII. Characterization of the Random SIR 

In this section we characterize the system performance as 
a function of the random quantity the signal to interference 
ratio (SIR) given by 

S{t) 



SIR(r) 



J(r)- 



(38) 



The SIR is a useful quantify for performance analysis in cel- 
lular systems because performance at the cell edge is usually 
interference limited; we consider additive noise and signal 
to interference plus noise ratio (SINR) in the simulations. 
Because simplified expressions including the exact distribution 
calculated in Section [rv] -Section[V]are only available in special 
cases, we use use the Gamma approximation described in 
Section |VI] To make the results concrete, we assume that 
S{r) = Sr where Sr ^ r[fcs,6's]. Note that in the distant- 
dependent path-loss and transmit power are absorbed into the 
9s term for ease of presentation. 

We use the SIR(r) to evaluate two metrics of performance: 
the success probability defined as P(SIR(r) > T) where T is 
some threshold and the ergodic achievable rate given by 



r(r) :=Ehi(l + SIR(r)). 



(39) 



The success probability is one minus the outage probability 
and gives a measure of diversity performance and severity of 
the fading. The achievable rate gives the average rate that can 
be achieved at r assuming that the interference is Gaussian. 

A. Success Probability 

We first provide the success probability without proceeding 
to a Gamma approximation for the out-of-cell interference and 
signal distribution given by S{r) = Sr where Sr ^ ^[ks, ds]- 
Using the results from | ,41J , the success probability at location 
r is given by 



P(SIR(r) > T) 



£{r)T 



' d^Ci{s) 



,_f(r)T 



(40) 



where in the case of heterogeneous interference, £/(s) is given 



by (16 1. In the special case of fc^ = 1 the success probability 
can be simplified as 



P(SIR(r) > T) = £,(£(r)r/(P,0,)). 



(41) 



Evidently, the closed-form expression ( |40| ) is involved and 
no useful insight can be obtained. This motivates the use 
of Gamma approximation for the aggregate interference as a 
means to provide simpler and useful expressions. 

We evaluate below the success probability for the case of 
Gamma distributed interference, effectively J{r) ~ T[ki,9.i\. 
In the case of homogeneous interference, the expression for hi 
is found in pTj i and for 9i in ( (32] i. In the case of heterogeneous 
interference, the expressions for ki and 6i are found in ( (35] l. 
The result is summarized in the following proposition. 

Proposition 15 (Success Prob. for Gamma Interference): 
The success probability when the signal distribution is 
S{r) — Sr where Sr ^ T[ks,9s] and the interference 
distribution is J{t) — Ir where ~ r[fci, 9i] is given by 



P(SIR(r) > T) = 



r(fc, + fc,) 
r(fc,) 



re. 



X 2F1 ki,ks + ki,l + h, 



TO, 



(42) 



where 2F1 is a regularized hypergeometric function. 

Proof: First we rewrite the probability to separate the sig- 
nal and interference terms as P(SIR(r) > T) = P (S*^ > TIr). 
Now conditioning on the interference and evaluating the 
expectation P(SIR(r) > T) = E{Fs{TIr)) = 



Tiks + ki) ( e. 



T{ks) 



TO, 



2F1 ki, kg + fci, 1 + ki,' 



T9, 



(43) 



where Fs{-) is the Gamma cumulative distribution function 
of Sr and the second equality follows from evaluating the 
expectation. ■ 
We can obtain a more exact expression in the case where 
each interferer is separately assumed to have a Gamma distri- 
bution. 

Proposition 16 (Heterogeneous Success Probability): The 
success probability when the signal distribution is S{r) ~ Sr 
where Sr ~ r[fc5,0s] and the interference distribution is 



J(r) = 4''^ 



J, 



(2) 



Jr'^'^^ where Jr™"* 



, , M is given by 

r{k, + p + n) f 



for m = 1 , 2 

n— ^ ' ^ 

X 2F1 [ n + p,ka + n + p\l + n + p\ — 



(44) 

where the parameters C, p, 9min are found using the results 
of Theorem |4l 

Proof: First we rewrite the probability to separate the 
signal and interference terms, and then conditioning on the 
interference, we evaluate the expectation, i.e. 

P(SIR(r) >T)^^(^r> tY^iIJ^^^ = E (^s{TY^l'r'>)^ . 

(45) 
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From the results of Theorem [4] the distribution of J2m 
can be expressed as 



or(P + n)CT" 



where 9^ 



3(m) 



_ Y^M , (m) ^ _ 

, and Cm is found by recursion from 

n+l 

(47) 



n+l 



with 7. = J:^^, k^'p-Hl-ermn/er'T- Now recognize 
that (|5|l is essentially a mixture of r[/ci + n, 6'inin] random 
variables with weights Cc„. Now we can use the results of 
( [T5| l to establish that 



p(siRM>r) = cf:c„£(%^t^(, 



n+p 



X 2-Fi (n + p,k3+n + p;l + n + p: 



Using this result, the average rate with homogeneous interfer- 
ence follows. 

Proposition 18 (Homogenous Average Rate): The average 
rate with homogeneous interference in the absence of noise 

is T(r) =Eln(l + SIR(r)) = 

C ^ c„ lrL{6„iin) + C* X! ^^'^^(P + '^) - i'i^i) - 

(50) 

where dmin Tnin{0s,9i\, p = + ki, C = 

{Q,ninlQs)^'(fi7mn/Qi)^\ and c„ is found by recursion in 

Proof: Follows from the application of Proposition [TT] to 
the case of the sum of two Gamma random variables r[fcs, 9s\ 
and r[fci,0i], and using the result in Lemma|2] ■ 
The average rate with heterogeneous interference can be 
computed along the same lines as Proposition [Ts] 

Proposition 19 (Heterogeneous Average Rate): The 
average rate with homogeneous interference in the absence 
of noise is 



00 



(48) 



While the result in Proposition [16] is more exact, it requires 
truncating the infinite series. In related work fyT\ we found 
that in most cases the approximation outperforms the truncated 
expression in ( [T7| ) unless a large number of terms are kept, 
and therefore conclude that the expression in Proposition 16 
is most convenient. 

The regularized hypergeometric function 2F1 (a, b, c, z) is 
available in many numerical packages or can be computed 
from unregularized hypergeometric function by recognizing 



r(r) = C)_Jc„ln(^?^:i) + C J] c„^(p(-) + n) 

n=0 n^O 

00 00 

D dn HelL) - ^ E + n) (51) 



n=0 



n=0 



Where := min{0„ \ \ . . . , '^}, p(-) = fc, + 



min 



:= min{0«,0f\...,0(^^)}, C 



kr' + --- + ki 

(t^mm/t^sj " llm=l(^miTi/t^i )' ' P^ - k^ H h fc, , 

^ = n^LiC^mL/^l'"^''"'. and Cm and are found by 
recursions like in (|6|. 

Proof: The first pair of terms result from application 



that 2F1 {a,b,c,z) = 2F^{a,b^,z) ^ p^^. ^^^^jj ^^j^^^^ of |z| it of_^ Proposition |17| to the term ElnS, 



r(c) 

is convenient to compute it using the truncated series definition 
[42 J , 1, 43 J while for larger values of |z| the asymptotic expres- 
sion can be employed | [44) . Depending on which numerical 
package is employed, it may be faster to compute the SIR 
using Monte Carlo techniques based on r[A:s, 6*3] and r[ki,9i]. 

B. Average Rate 

The average rate is useful for computing average rates at the 
cell edge and the area spectral efficiency, where the average 
rate at r is integrated over the radius of the cell. Our calcula- 
tions are be based on the observation that E ln(l + SIR(r)) = 
Eln(S'r + Ir) - Eln(/,.). From Lemma |2] Eln(/^) has a 
computable solution but Eln(5,- + Ir) requires dealing with a 
sum of Gamma random variables. To solve this problem, we 
can use the result in (|5]l to find an expression for the expected 
log of the sum of Gamma random variables. The resulting 
expression is summarized in the following proposition. 

Proposition 17: Suppose that are independent 

r[fcm,^'m] distributed random variables. Then 



, while the second term results from appli- 



(M) 



cation of Proposition [it] to the term E In l^'' + /r'^' • • • /j 



If the number of terms is too large for practical computation, 
an alternative is to approximate the sum of the signal plus 
interference term as another Gamma random variable using 
the moment matching approach. 

Proposition 20 (Homogeneous Average Rate - Gamma): 
The average achievable rate without noise with homogeneous 
interference is approximately 



T(r) K, ip 



In 



ks9s + kz9z 



k,9l + kz9l 

- i^{ki) - \n{9,) (52) 

Proof: To find Eln(S',. + Ir) we apply Lemma |5] and 
approximate the sum of Gamma random variables with another 
Gamma distribution with the same mean and variance. Note 
that 



00 n— 1 



E S'r + Jr = ks9s + fci^j 

var Sr + Ir = ks9^ + k^9f. 



(53) 
(54) 



E In Xm = ^(p) + ln(0m,„) + cY,CnJ2^f 

\m=l J n=0 m=0 ^ * 

(49) 

Proof: See fTP, Proposition 5]. ■ 



Then a Gamma random variable Z ^ ^[k^, ^z] with the same 
mean and variance has 

{k,9s + kz9z)^ ^ _ ks9l + kz9l 



ks9^, + kz9l 



(55) 
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Now we can approximate the term 



In 



(56) 



The last part of the expression for E In 1^ follows from Lemma 

m ■ 

The result can be extended to the case of heterogeneous 
interference by modifying the computation of the mean and 
variance terms. 



VIII. Simulations 

In this section we consider a system model of the form in 
Fig. [2] with a single antenna base station and a Poisson field 
of base station interferers. We evaluate performance in a cell 
of radius Rc = 300m with a guard region of Rg — 150m. 
The base station locations are modeled as a PPP with A ~ 
2/\/27i?c = 4.3 10^^ transmitters per which corresponds 
to the same density as a hexagonal lattice with cells of radius 
Rc- The base station has 40 transmit power with a single 
antenna. Independent identically distributed complex Gaus- 
sian fading is assumed (corresponds to Rayleigh distributed 
channel amplitudes). The path-loss model of ( [T3] l is used with 
n — 3.76, do — 35m, and C = 33.88. Log-normal shadowing 
is assumed with 6dB variance. Monte Carlo simulations are 
performed by simulating base station locations over a square 
of dimension 15Rc x 15i?c- The mutual information and the 
outage probability are estimated from their sample averages 
over 200 small scale fading realizations, and 500 different 
interferer positions. A probability of success threshold of 
T = 31 (corresponding to 5dB) is used. Performance is 
evaluated as a function of the distance from the center of the 
fixed cell. 

First we validate two of the key assumptions used in most 
of the results in this paper: the use of a guard region and 
the approximation that the interference power for a user at 
radius r from the cell center is upper bounded by considering 
interferers outside a ball of radius Rc + Rg — r. In Fig. |4] 
we compare the performance between the proposed hybrid 
approach and a hexagonal layout. We consider three different 
guard regions for the hybrid approach: Rg = 0, Rg = Rc, and 
Rg — 2Rc/3. We compute the received interference power 
for a user at radius r from the cell center by considering that 
the interference lies outside a ball of radius Rc + Rg — r as 
illustrated in Fig. |2] We compare with a hexagonal layout of 
cells with either one tier of interferers (six total) or two tiers of 
interferers (eighteen total). We plot the ergodic rate in Fig. |4] 
and the probability of success in Fig. [5] Our approach with a 
guard region of Rg — 2Rc/3 matches the multi-tier hexagonal 
layout both in terms of ergodic rate and probability of success. 
Consequently our numerical results justify our claim that a 
Poisson interference field can approximate an interference field 
from base stations located on a hexagonal grid if a suitably 
chosen guard region is employed. 

Now we consider the Gamma approximation for the inter- 
ference field. We claim that the Gamma distribution is a good 
approximation of the distribution strength of the interference 
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R Guard 
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Hexagonal layout: 2 tiers of interferers 
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Fig. 4. Illustrates the ergodic rate as a function of radius from the cell 
center. The flat part at the beginning is due to the non-singular path-loss 
model that is constant until the reference distance of do- The case of no 
guard interval provides a lower bound on performance while a guard interval 
of Rc provides an upper bound on performance. Performance without the 
inteii'erence field approximation is good near the cell center but deviates 
further away from the base station. The hexagonal layout with one and two 
tiers of interference almost exactly overlaps the performance with a guard of 
2Rc/3. This justifies both the symmetric interference field approximation and 
the use of a guard interval Rg = 2Rc/3. Note that in this case we compare 
the Poisson distributed interferers with the hexagonal layout of interferers to 
justify our claim, there is no Gamma approximation of the interference. 
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Fig. 5. Probability of success curves derived from the Kaplan-Meier estimate 
of the cumulative distribution function at a distance of r = SOOm. Due 
to the number of points this curve is best viewed in color. Similar results 
are observed as for the average rate. The case of no guard interval provides 
a lower bound on performance while a guard interval of Rc provides an 
upper bound on performance. Performance without the interference field 
approximation is good near the cell center but deviates further away from 
the base station. The hexagonal layout with one and two tiers of interference 
almost exactly overlap the peiformance with a guard of 2Rc/3. This further 
justifies both the symmetric interference field approximation and the use of a 
guard interval Rg = 2R.c/3. Note that in this case we compare the Poisson 
distributed interferers with the hexagonal layout of interferers, there is no 
Gamma approximation of the interference. 



power. As an error metric we use the percent error in the 
expected log as given by 



|ElnX„ 



error 



Eln/(r) 



|Eln/(r)| 
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The motivations for comparing ElnX are that it shows up 
in the calculation of the average rate and that signal and 
interference power are usually measured on a log scale. Note 
that we match moments on the linear scale (not the log scale) 
thus the comparison seems to be fair. The comparison is 
described in detail in Fig. [6] Based on our comparison we 
conclude that with a suitably chosen guard region of 2Rc/3, 
the gamma approximation gives an error in the log that is 
acceptable, justifying our use of the Gamma approximation of 
the interference field. 




- base stati 

- base stati 

- base stati 

- base stati 

- base stati 
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ons witti 10m radius 
ons witti 30m radius 
ons witti 60m radius 
ons witti 90m radius 
ons witti 120m radius 
ions witti 150m radius 
ions witti 300m radius 
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Fig. 6. In this simulation, tlie intensity of tlie interference field is varied 
according to A = 2/ \/27 by changing the size of Ri . There is no guard 
region. It can be seen that the approximation is better for higher intensity 
interference. For our simulations with Rc = 300m we see that the percent 
en'or at the location of the guard region 2Rc/3 = 200m is less than 3%. 



probability, similar results are observed. There is some error 
with the homogeneous interference at the cell edge, the fit 
with homogeneous and DAS interference is good, and the 
inclusion of femtocells reduces the accuracy slightly. While 
there is approximation error, we still believe that the fit is good 
enough to justify the viability of the Gamma distribution. 

Heterogeneous Network 
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Fig. 7. Comparison of the ergodic rates with and without the Gamma 
approximation with a tier of low-power femto nodes that create cross-tier 
interference. We see that the Gamma approximation provides low en'or except 
towards the cell edge for the case of homogeneous interference. With cross-tier 
interference, the Gamma approximation provides a constant offset of the entire 
simulation curve. The offset is constant because the cross-tier interference is 
constant. The en'or though is consistent, meaning that it is still possible to 
make relative comparisons between the Gamma approximation curves. 



Now we compare the Gamma approximation in terms of 
the average rate and the success probability in a case with 
heterogeneous interference. We consider a setup with a single 
antenna base station and three sources of interference. One 
source of interference is from other base stations, which 
has power and density as we already described. The second 
source of interference is a multiuser MIMO distributed antenna 
system (DAS) with four single antenna transmission points, 
a single active user per cell, and the same total power The 
DAS system is also modeled using a PPP but with density 
4A — 1.7 10~^. The distribution of the effective channels 
are derived from |37|. Note that we assume essentially that 
every cell also has a DAS system, in reality some cells would 
have DAS and some would not, which can be included by 
modifying the density. The third source of interference is 
from cross-tier interference from a single tier of femtocells. 
The femto interferers have a 2QdBm transmit power and a 
guard region of 5m. Since this is cross-tier interference, the 
femtocells may be deployed in the cell and thus they create 
a constant level of interference power corresponding to a 
minimum distance of 5m. An indoor-outdoor penetration loss 
of 20dB is assumed. The results of the comparison in this 
section are papered in Fig.[7]for average capacity and in Fig. [8] 
for outage. In terms of average rate, the main conclusion is 
that the Gamma approximation provides good performance for 
heterogeneous interference. There is a more significant error 
at the cell edge for only a single interferes In terms of success 



Heterogeneous Network 




Radius 

Fig. 8. Comparison of the success probability with and without the 
Gamma approximation for a target SIR of 15dB. We see that the Gamma 
approximation overall provides reasonable performance except right at the cell 
edge. At small SIR it tends to underestimate the outage. There is about a 15% 
error with small cell interference versus having homogeneous interference. 
The outage curves were computed using Monte Carlo techniques. 



IX. Conclusions 

In this paper, we proposed a hybrid model for determining 
the impact of interference in cellular systems. The key idea is 
to develop a model for the composite interference distribution 
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outside a fixed cell, as a function of the user position from 
the cell center. From a numerical perspective, our approach 
simplifies the simulation study of cellular systems by replacing 
the sum interference term with an equivalent interference 
random variable. Then functions like average rate and success 
probability can be computed through numerical integration 
rather than Monte Carlo simulation. Unfortunately, the numer- 
ical integrals may still require a fair amount of computational 
power to compute. Consequently, we proposed to approximate 
the interference distribution by moment matching with the 
Gamma distribution. We showed how the Gamma distribution 
could be used to obtain relatively simple expressions for the 
success probability and ergodic rate, which simplify further 
when the sum interference power is approximated directly as 
a Gamma random variable. Simulations showed that with an 
appropriately chosen size of the guard radius, a reasonable fit 
between the case of PPP distributed co-channel interference 
and hexagonal interference was achieved. Further work is 
needed to better characterize the approximations, e.g. by 
developing expressions for bounding the error terms. 
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